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0> . Abstract 



^ ■ The dynamical content of pNN and pNA vertex functions is studied with a mesonic 

model. A set of coupled integral equations satished by these vertex functions were 
solved self-consistently. These solutions indicate that the dominant mesonic content 
arises from di-pion dynamics. With the experimentally determined pion-baryon-baryon 
coupling constants and ranges as input, the model predicts a gpNN that agrees with the 
meson-exchange-potential results. On the other hand, it predicts a smaller fpNA and 
much softer form factors. Implications of the findings on the use of phenomenological 
coupling constants in nuclear reaction studies are discussed. 
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^ ■ PACS index: 14.20.-c,Gk; 21.30. +y; 25.80.-e 



Meson-exchange theories have been very successful in describing nuclear phenomena 
in the region of low and intermediate momentum transfers, where practical methods for 
solving non-perturbative QCD are yet to be fully developed. Within the framework of the 
meson-exchange theory, understanding the dynamics associated with the pNN and pNA 
vertices is a subject of great interest. It is well known that p-exchange gives rise to a nucleon- 
nucleon (NN) tensor force which has a sign opposite to that arising from the exchange of 
TT. The partial cancellation between these tensor forces makes it indispensable to include p 
in the hadronic description of nuclei and nuclear matter. The important role of p in nuclear 
structure has further motivated the study of its role in nuclear reactions. Studying the p 
is, however, complicated by the fact that, unlike vrA^ — > A, the pN -^ N and pN -^ A 
processes do not occur in free space. Consequently, information on these reactions were 
mainly obtained from fitting NN phase shifts with the use of meson-exchange potentials 
(MEP). It is, therefore, quite likely that the MEP parameters may have to be modified for 
studying processes other than A^A^ scattering. Indeed, analyses of deep inelastic scattering 
data indicate that the range of a pion monopole form factor cannot exceed 650 MeV. f!]] |@| 
This upper limit is nearly a factor of two smaller than the MEP ranges. Smaller ranges 
have also been predicted as a necessary outcome of the boson nature of the pion. A recent 
study of the p{p,n)A^~^ reaction has shown that the inclusion of p-exchange mechanisms 
actually worsens the fit.|^ This disagreement raises the interesting question as to what are 
the appropriate p-meson parameters to be used in meson production. These developments 
have motivated us to carry out a detailed analysis of the pNN and pA^A vertices by means 
of a dynamical model that does not treat the rho-baryon-baryon {pBB') processes as contact 
interactions. The main goal of our study is to identify important dynamical contents of the 
phenomenological coupling constants. 

Our model for the pNN vertex is illustrated in fig.l, where the wavy, dashed, thin, 
and thick lines denote, respectively, the p, vr, A^, and A. The initial and final nucleons are 
labelled a and c, while the intermediate baryons are labelled b and b'. Figure 1 corresponds 
to the equation 

Fpa;c = Fp.^TnTG)j._^' G\^,J FT,-^a;bG i^ -^7r2fe;c + Fa-^nbGi, Fph-fi^Gf,, F-jrh'^c = Bpa-c + Fpa;c ^ (1) 



where Fpa-^^ denotes the pa —>■ c vertex function. A similar notation apphes to the other vertex 
functions, with the indices b, V representing either an A^ or a A. In Eq.([I|) Bpa-^ denotes the 
Born amphtude given by the two triangle diagrams in fig.l, while F' corresponds to pion 
corrections of the ph -^ h' vertices, as shown by the four vr-loop diagrams. Since the solutions 
of the model, Fpi-j {i,j = a,b,b'c), appear on both sides of Eq.(0), we must solve a set of 
four coupled equations of the form of Eq.([l|) with the baryon indices ac equal to A^A^, A^A, 
AA^, and AA, respectively. For succinctness, only the equation corresponding to ac = NN 
is illustrated in fig.l. We have solved the coupled equations by iteration, and a convergence 
was obtained after five iterations. In principle, one should also include the p-loop corrections 
of the vertices. Because Mp » M,r, these corrections will be much less important than the 
vr-loop corrections which, as we shall see, are quite small. They were thus neglected in our 
calculations. 

It is worth noting that diagrams similar to ours have recently been considered in ref.0. 
We emphasize, however, that there exist important differences between the analysis in ref. |Q 
and ours. In ref. [^ the leading contributor was assumed to be the phenomenological three- 
branch pNN vertex while the triangle and loop diagrams were considered as corrections 
and treated perturbatively. Moreover, phenomenological MEP pBB' coupling constants 
were used as inputs and no coupled equations were solved. As such, it is equivalent to 
evaluating only first-order corrections to the MEP pNN coupling constant. Consequently, 
the calculation does not address the djTiamics leading to the phenomenological coupling 
constants. In contrast, the analysis presented in this paper considers the triangle diagrams 
as the driving terms and the loop diagrams as corrections. By solving the coupled equations, 
the pBB' coupling constants and ranges are obtained as the self-consistent solutions of the 
model considered. Hence, these solutions contain valuable information on the dynamical 
content of the coupling constants. 

We used the following partial-wave decomposition for the pBB' vertex functions: 

SmsLrriL 

X < Sms,LmL \ {SL)JcVc > Yl^^{p)J^(^pp,pa,Pc) , (2) 



with the radial vertex function JF parametrized by 

^{Pp,Pa,Pc) = G{L)pa;c^==P^VL{Apa;c,p) ■ (3) 

The Pi {i = p, a, c) denote the four-momenta of the i th particle, and pi is its spatial 
part; p denotes the magnitude of the relative momentum between p and a, and Wc is the 
invariant mass of c. The two Clebsch-Gordan coefficients specify, respectively, the isospin 
and angular momentum coupling schemes with I and J denoting the isospin and spin of the 
particles, and t, v their z-components. The vertex function JF is a scalar and depends on two 
independent four-momenta which can be chosen as Pp + Pa and Pp — Pa- It is convenient to 
work in the p — a cm. system where the external four-momenta have the simple expressions 
Pp = (p°,p). Pa = (a°, — p), and pc = {wc,0). Clearly, Wc = a° + p°. From the two 
independent four-momenta one can form three independent invariant scalars. We will put 



the baryon a on the mass shell so that a° = Ea{p) = J ml + p^. Consequently, JF will 
depend only on two independent variables which we choose as Wc{= w) and p{=\ p |). 
Because the parametrization used in Eq.(^ is not the most general one, we should expect 
Q and A to depend on w. This w-dependence has important physical consequences and will 
be discussed later. In the literature, vertex functions have often been made to depend only 
on one variable. We emphasize that the one-variable dependence is exact only when two of 
the three particles are on their mass shells. 

Parity conservation limits the values of orbital angular momentum in Eq. (^ to L = 1 
for pN -^ N and pN — * A, and to L = 1 and 3 for pA -^ N and pA -^ A. As the Lagrangian 
models in the literature consider only L = 1, we shall solve the coupled equations in the 
p-wave channel and omit, henceforth, the index L oiQ and v. Parity conservation also limits 
the p ^ irn, ttN -^ N, nN — ^ A, and vrA -^ N processes to relative p-wave interactions 
alone. On the other hand, ttA -^ A can have both p- and f-wave interactions. Again, 
only the p-wave interactions will be retained in the calculations in order to make a close 
connection to the Lagrangian models. The ptttt vertex function is parametrized as 



Fp-^^{Pp,PuP2) =Y.'^iH—p==KAp^^^^)^im{k)^2uj^{pi)2uj^{p2) , (4) 

where Ci =< Iti, 1^2 | (ll)Wp >,h = kr/{l + k'^r'^Y with r = l/Ap^^r, and fc =| k | is the 



TTTT relative momentum. The TTiab vertex function is parametrized by 

Fma;b{PuPa]Pb) = J2 C iC j^^^^^h{Ka;h, K)Y*^^{k) , (5) 



where Cj =< lata, Iti | {Jal)Ibtb >, Cj =< Ja^a, l"^i | {Ja^)Jb^b >, 1^ IS the Tia relative 
momentum, mi is the z-component of the corresponding orbital angular momentum, and h 
is a dipole form factor defined by h = (k/Mtt) [(A^ — Mtj'^)/{K^ + k^)]^. The vertex functions 
F-K2b;c, Fpb-b', Fa-nb, and F.„b/-c are parametrized in terms of the corresponding momenta in 
a form similar to Eq.(|^). The triangle as well as the loop diagrams depend on one four- 
momentum integration variable q = (g°, q). By closing the contour of q^ along a semi-circle 
in either the upper or the lower half-plane, we can carry out the g° integration analytically. 
Upon introducing Eqs.(0)-(|^) into Eq.(P and projecting out the (L = 1) angular momentum 
dependence, we obtain a set of four coupled equations for v. 

One should note that the value of a coupling constant depends on the parametrization 
convention of a theory. It is, therefore, useful to relate the ^'s of the present partial-wave 
formalism to the coupling constants of other works. To this end, we have used the S- 
matrix convention of ref . [B] . In relation to the Lagrangian model and the parametrization 



employed in ref.|0], we obtain G-^-nn = yi3/2iT^)fj,NN and GttNA = /ttAta/vStt^- For the 
ttAA vertex, ^^aa = (5/2)/;^AA/v^67r2 where /^^A = {M^/2MA)gnAA and ^f^AA is the 
axial- vector coupling constant defined in ref. [Q. With respect to the parametrization in 



ref.0, gp^^/^2Mp = Qp^^ and ^^aa = (15/4)/^AA/v^67r2- Hence, /^aa = (2/3)/4aa- In 
addition, f-^NA and gp.„.„ can be directly related to the experimental widths of A and p. 

The QpBB' can also be readily related to the /'s and f^'s defined in the Lagrangian 
formalism. Upon expressing the matrix elements < B' \ TipBB' I pB > of the Lagrangian 
model in a partial- wave representation and equating them to Eq.(^, we obtain at p^ = 

iGv)pNN 4(^2 - l)y'37r2MpMiv = {fu)pNN + {gu)pNN = (1 + K,p){9u)pNN , (6) 

{gv)pNA)Jl{V5 + l)^37rm^p/M^ = (/n),^A, (7) 

where Qv has the dimension of inverse mass, and Hp = {f/g)pNN- The u denotes the specific 
form factor employed in the Lagrangian model, which varies from one published work to 



another. Using these two equations, one can make a straightforward comparison between 
the di-pion and the MEP results. Our calculation gives the sum (/ + g)pNN- A knowledge 
of Up from the analysis of the nucleon electromagnetic form factor can be used to determine 
/ and g separately. 

For the calculation, the values Qp^^T, = 0.6684M~^/^, Ap^rTr = 2.336 fm~^, fj^NN = 1, 
and A^ATjv = ^nNA = 950 MeV(= A(2)) were used. The pirn parameters fit the vrvr phase 
shifts. The value of the dipole range A(2) is taken from ref. 0. It corresponds to a monopole 
range A(i) = 0.644 x 950 = 612 MeV, which is nearly half of the MEP pion monopole range 
of 1.2 GeV. [0[|lO| For the other parameters, we used the following sets: (a) fnNA = 1-69 



(ref.0]), /^AA = 0.8 (ref.i), A^aa = A(2); (b) Una = 2.19, /,aa = 1.03 (or /^^^ = 1.55), 
AvrAA = A(2); and (c) fnNA = 2.19, fnAA = 1-03, A^^aa = 1-86 GeV corresponding to the 
extreme situation of having an equivalent monopole range A(i) ^^aa = 1-2 GeV. The criteria 
leading to the above choices are as follows. In set (a) the ratio fnAA/ fnNA = V^/3 satisfies a 
SU(6) quark model result. |TT| Using the experimental Fa = 115 MeV to calculate fnNA and 



using gl/^^/4:7i = 60 (given by the isobar analyses [|]|]T2| of the iiN — > vrvrA^ data), we obtain 
set (b). As these isobar analyses do not give A,rAA, we used A^tAA = 950 MeV in set (b) and 
1.86 GeV in set (c) in order to see the effects of this range parameter. We have also evaluated 
contributions by intermediate A^*(1440) state in the triangle diagrams, using fnNN* = 0.467 
and fnAN* = 1.63 calculated from the median values of experimental partial widths TN*{nN) 
and TN*{nA)- The full, coupled-equation results are presented in Table 1. The values within 
parentheses are due to the inclusion of the A^* contribution. The (/ + g)pNN and fpNA have 
been evaluated, respectively, at Wc = 939 and 1232 MeV; they are found to be higher than 
the corresponding Born results by 5% and 14%. For the convenience of making comparisons, 
we have listed the values of / and g, as obtained from Eqs.(P) and (^. Further, gpNN has 
been calculated from the solutions (/ + g)pNN with Kp = f /g = 3.7 given by an analysis of 
the electromagnetic form factor of the nucleon. 

The MEP values of gpNN vary from 1.28 to 2.27 (i.e. g^/An = 0.13 - 0.41). Those 
for fpNA vary from 4.91 to 7.81 (or fpNA/'^'^ — 1-92 — 4.86). These lower and upper limits 
are, respectively, the results of refs.[^] and ^. An inspection of Table 1 shows that with the 
inclusion of A^*, sets (b) and (c) lead to a gpNN that is very close to the MEP values. In 
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fact, a good agreement has been obtained when we used f-rrNN* given by the experimental 
upper bound of TN*{TrN)- While parameters of set (a) give a real part of fpNA that is smaller 
than the lower limit of the MEP values by a factor of ~ 2.5, the use of fnNA and /^tAA 
derived directly from the data (set (b)) brings the difference down to within 40%. Use of a 
large A^aa (set (c)) makes the real part of f^NA agree with the MEP values. However, in 
view of the deep inelastic scattering data on the upper limit of A^jvjv and A^^ta, the large 
AttAA of set (c) may be questionable. Consequently, we consider the results due to set (b) as 
being more realistic. Our results also indicate that the contribution by N* is small. At this 
point, a comment on the complex nature of fpNA is in order. As pointed out after Eq.(||), 
the coupling constants can be energy- dependent. An important feature of the di-pion model 
is that fnNA becomes complex-valued at w > Mj, + Mat = Wth- The fpNA of Table 1 
were calculated ai w = M^ > Wth, where the vrA^ channel is open, thereby, giving rise to 
Im[/pAr:A] < 0. The energy dependence of the real and imaginary parts of the pNA form 
factor at p=0 are shown, respectively, as the solid and dashed curves in fig. 2, where we have 
defined H{w,p) = {Qv){w,p). 

Table 1: The coupling constants and ranges (in MeV) given by the di-pion model 

for the pNN vertex a.t w = Mn and for the pNA vertex at w = Ma- 



Set 


(/ + 9)pNN 


QpNN 


ApATAT 


fpNA 


ApATA 


a 


4.39(4.90) 


0.934(1.04) 


435 


1.89-2.09i(2.08-2.09i) 


495 


b 


4.82(5.32) 


1.025(1.13) 


430 


3.50-2.7H(3.69-2.7H) 


490 


c 


4.86(5.36) 


1.033(1.14) 


430 


5.25 - 2.70i(5.36 - 2.70i) 


490 



Table 1 indicates that the ranges of the calculated form factors are in the region of 
450 MeV, much smaller than the MEP values of 1.2 GeV. |^[|T^ Here, we are considering the 
form factors as functions of p^ and define A as the momentum at which the magnitude of 
form factor is half of its value at p^ = 0. Similar ranges have been obtained when we used 
the four-momentum transfer p^ as the variable. We can, in fact, easily show that a smaller 
range is a consequence of the composite nature of a vertex. 

We have further examined the four-pion content of the vertex functions by adding 
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the p -^ UJTT mechanism into the triangle diagrams. Using parameters of set (b) together 
with g^p^M^=9.71 (ref.]TB[), 5f^ArAr=11.54 (ref.0]), and assuming the relation QujAA/gujNN = 
Q-kAaI Q-kNN in our calculations, we have found that contributions from four-pion dynamics 
are negligible because M^^ ^ M.,^. Although diagrams having more than three subvertices 
could in principle also contribute, they tend to be unimportant because more are the particle 
propagators less is the interaction probability. It is, however, worth recalling that while both 
the L = 1 and 3 interactions can contribute to the ttAA process, only the L = 1 case is 
considered in this and other works. Although inclusion of the L = 3 interaction can be easily 
implemented in the present formalism, experimental information on this f-wave coupling is 
sparse and uncertain. If one can establish experimentally that the f-wave coupling constant 
is not too small, then it could make a sizable contribution because the f-wave vertex function 
is oc q^ and q is an integration variable that extends to very large values. 

In summary, the di-pion model predicts QpNN that agrees with the MEP values, but 
the predicted fpNA is about 30% smaller. In view of the fact that the QpNN and fpNA of 
this work are based on a microscopic model while those of MEP were correlated fitting 
parameters, we regard the results given by the two approaches as being compatible. The 
inclusion of four-pion dynamics via the ujti doorway state does not alter the results, indicating 
that the di-pion mechanism does represent the leading mesonic contribution. However, 
experiments capable of providing information on f-wave vrAA couplings can help determine 
if mesonic dynamics alone is sufficient to account for all the strength of the phenomenological 
pA^A coupling constants given by MEP. The ranges of the calculated form factors are in the 
450 MeV region, much smaller than the corresponding MEP values, but in line with the 
small ttNN and vrA^A ranges. |[I|] [|2[ As mentioned earlier, a small range also reflects the 
composite nature of a vertex function. Smaller ranges will give calculated meson-exchange 
contributions that are different from those given by large ranges and, hence, may lead to 
new understandings of the data. Our analysis has shown that in the medium-energy regime, 
there is no compelling need for introducing non-hadronic dynamics. We stress, however, that 
MEP theory will become impractical at very high energies because the exchange of many 
heavy mesons has to be included in order to produce the correct energy dependence of the 
A^A^ cross section. Furthermore, the theory will be deficient because the quark substructure 



of the hadrons will start to manifest. We also stress that while it is a good approximation to 
employ real-valued vertex functions in analyzing A^A^ scattering below the pion production 
threshold, |T^ the situation is different in meson production experiments where the energy 
Wth can be surpassed. Thus, using real-valued vertex functions is questionable. We beheve 
that a nonvanishing imaginary part of the form factor can give new interference effects 
in nuclear reaction calculations. This aspect of the meson-exchange dynamics merits a 
systematic investigation in the future. 
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Figure 1: Feynman diagrams for pBB' vertex functions. 
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Figure 2: Re[if (w, 0)] (solid curve) and —Im.[H{w, 0)] (dashed curve) of the Born pNA vertex 
function versus {w — M]\j)/Mt^ obtained with parameters of set (b). 
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